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LONG-TERM  GOALS 

Radio  tracking,  communications  and  guidance  at  sea  rely  on  interpreting  of  radar  signals.  This  has 
motivated  extensive  studies  of  the  physical  factors  governing  the  pattern  of  signal  propagation  over  the 
ocean.  Although  over  the  last  several  decades  the  models  of  electromagnetic  propagation  have  been 
substantially  advanced  (Kukushkin  (2004);  Levy  (2000)),  measurements  and  modeling  results  still 
show  differences  that  can  distort  signal  interpretation.  In  particular,  signal  intensity  at  the  receiver  is 
often  over-predicted  by  models  (Barrios  and  Patterson  (2002)),  which  suggests  that  mechanisms  of 
signal  energy  loss  are  still  incompletely  understood  and/or  are  not  correctly  built  into  the  propagation 
models.  The  long-term  goal  of  this  work  is  to  identify  the  source  of  discrepancies  between 
observations  and  propagation  model  predictions  and  update  the  models  so  that  the  discrepancies  are 
reduced. 

OBJECTIVES 

An  essential  element  in  models  of  over-the-ocean  propagation  is  the  scattering  from  the  wavy  water 
surface,  a  problem  we  specifically  focus  on  in  this  work.  Although  the  understanding  of  rough  surface 
scattering  has  been  gradually  expanded  (Elfouhaily  and  Guerin  (2004);  Voronovich  (1999))  and  the 
exact  solution  of  the  scattering  problem  can  be  obtained  numerically,  computational  efficiency  requires 
that  we  use  a  simplified  description  of  the  scattering.  For  that  purpose  one  specific  scattering  model 
(Miller,  Brown  and  Vegh  (1984))  has  been  considered  to  be  more  realistic  (Levy  (2000))  when 
compared  to  earlier  models  (Ament  (1953);  Eckart  (1953b))  and  has  been  adopted  in  propagation 
models  employed  by  the  US  Navy.  However,  an  analysis  of  that  model  determines  that  it  is  inconsitent 
with  preexisting  both  theoretical  predictions  and  observational  results  regarding  the  statistics  of  the 
ocean  surface. 
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Figure  1.  Probability  density  distributions  of  surface  elevation  (upper  left),  of  wave  amplitude  (lower 
left),  and  of  wave  phase  (upper  right.  The  experimental  data  were  obtained  over  1  hour  of  surface 
elevation  measurements  during  the  Rough  Evaporation  Duct  experiment  The  vertical  lines  in  the 
upper  left  plot  mark  the  waves '  standard  deviation.  The  filled  circles  indicate  experimental  data,  the 
continuous  line  corresponds  to  the  case  of  Gaussian  distribution  for  the  surface  elevations,  and 
dash-dot  line  shows  the  probability  distributions  as  assumed  in  the  Miller-Brown-Vegh  model.  The 
lower  right  plot  presents  the  magnitude  of  the  scattered  field  according  to  Ament  (1953),  for  a 
Gaussian  random  surface  (continuous  line)  and  according  to  Miller-Brown-Vegh,  dash-dot  line. 


The  model  of  Miller,  Brown  and  Vegh  (1984)  is  based  on  the  Kirchhoff  approximation  (i.e.  assuming 
the  curvature  of  the  rough  surface  to  be  much  greater  than  the  wavelength  of  the  electromagnetic 
signal).  The  model  also  prescribes  some  statistical  properties  of  the  sea  surface,  as  summarized  in 
Figure  1.  Specifically,  the  model  of  Miller,  Brown  and  Vegh  (1984)  considers  the  sea  surface 
elevation  as  described  by 


=  ,4  sin 


where  A  is  a  Gaussian  random  variable,  A  is  a  characteristic  wavelength,  and  x  is  unifonnly 
distributed  in  the  interval  (-A/4,  A/4) . 
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Given  that,  surface  elevation  7  is  a  random  variable  with  a  distribution 
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where  K0(  )  is  a  modified  Bessel  function  of  the  second  kind.  Within  the  Kirchhoff  approximation  the 

scattered  field  intensity  is  determined  by  the  characteristic  function  corresponding  to  the  distribution 
p(rj) ,  which  in  this  case  becomes 

R  =  exp(-2a,y  )/0(2cr^2). 


Figure  1  shows  the  probablility  density  distributions  of  surface  elevations,  wave  amplitudes,  and  wave 
phases  as  assumed  in  the  model  of  Miller,  Brown  and  Vegh  (1984)  and  for  the  case  of  Gaussian  sea 
surface,  as  in  Ament  (1953)  and  Eckart  (1953b).  Sea  surface  statistics  obtained  from  measurements 
during  the  Rough  Evapration  Duct  experiment  is  shown  in  Figure  1  as  well.  Clearly,  sea  surface 
statistics  proposed  by  Miller,  Brown  and  Vegh  (1984)  is  not  supported  by  the  observations.  The 
narrow  peak  in  the  distribution  of  elevations  appears  unphysical  and  has  never  been  reported  from 
measurements.  The  amplitude  must  remain  non-negative  to  be  consistent  with  its  definition,  thus  a 
Gaussian  distribution  cannot  be  valid  for  it.  The  phase  distribution  in  the  model  of  Miller,  Brown  and 
Vegh  (1984)  does  not  cover  the  phase’s  range  of  variation.  The  surface  statistics  inconsitencies  make 
the  model  of  Miller,  Brown  and  Vegh  (1984)  to  predict  enhanced  intensity  of  the  scattered  signal  and 
therefore  could  be,  at  least  partially,  responsible  for  the  tendency  of  propagation  models  to  overpredict 
the  signal  intensity  at  the  receiver. 

In  this  context,  our  objective  was  to  examine  the  model  of  Miller-Brown- Vegh  and  propose  a  revision 
that  is  free  of  that  model’s  deficiencies. 


APPROACH 


For  a  Gaussian  random  surface  r/(x,y )  ((7)  =  0  ),  all  the  central  moments  jUn  =0]")  can  be  expressed 
through  the  surface’s  standard  deviation  <j2  =  /u2  =  (r/2\  as 


(2») 


_2  n 


il2„  ~ WTT °"n  and  p2n+l=0. 
2  nl 


Specifically,  //4  =  3cr4 ,  ju6  =  15crb ,  //8  =  105cr5 ,  juw  =  945cr10 ,  etc 


We  used  data  collected  during  the  two  weeks  of  the  Rough  Evaporation  duct  experiment  to  obtain  the 
moments  of  the  sea  surface.  Figure  2  presents  the  probability  density  function  for  the  moments  of 
orders  between  2  and  16.  It  was  observed  that  the  moments  generaly  cluster  about  their  Gaussian 
values  although  do  not  follow  them  exactly,  i.e.  the  surface  statistics  shows  a  detectable  deviation  from 
Gaussianity. 
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The  analysis  of  cumulants  is  a  well  established  statistical  approach  that  has  not  been  applied  to  the 
problem  of  sea  surface  scattering.  We  chose  the  cumulant  analysis  as  an  adequate  tool  for  describing 
the  statistics  of  non-Gaussian  processes.  Such  approach  should  let  us  produce  an  approximation  for  the 
scattering  cross-section  of  the  sea  surface,  based  on  surface’s  statististical  moments,  which  are 
measurable.  Such  approximation  will  be  free  of  the  inconsistencies  inherent  to  the  model  of  Miller, 


Brown  and  Vegh  (1984).  where  the  cumulants  Kn  are  defined  by  Kn  =  A'("i(0)  =  i 
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Figure  2.  Experimental  probability  density  distributions  of  the  sea  surface  statistical  moments 
ju2,/u3  /  cr3,...,//l(i  /  cr16.  The  dashed  vertical  line  marks  the  value  of  these  moments 

for  a  Gaussian  variable,  Le.  (2n)\l{l" n>). 
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They  can  be  expressed  through  the  measurable  moments  as 


k2  =  Mi 

K3  =  Ms 

'U  =  Ma  ~  3/T 

*5  =  Ms  -^Mith 

«6  =  Mb  - 1 5 MaMi  - 1  °/^3  +  30//23 

^7  =  /T  -  2 1//5//2  -  3  5//4//3  +  2 1 0//,//; 


The  explicit  form  for  /rn  was  obtained  up  to  the  order  of  16.  For  a  strictly  Gaussian  random  variable, 

the  cumulants  of  order  greater  than  2  should  be  zero.  Figure  3  presents  the  distribution  density  of  the 
sea  surface  cumulants  calculated  from  Rough  Evaporation  Duct  measurements.  The  sea  surface’s 
deviation  from  Gaussianity  is  also  apparent  from  these  distributions. 

The  field  intensity  as  a  result  of  the  incoherent  scattering  from  the  sea  surface  is  determined  by  the 
joint  characteristic  function 


_  00  n  ## 

(p{qx ,q2)  =  \ ei(q'n' +q 272  )p{q] , / j2 )d q,d 72  =  Z  Z  ,,,  ''"Z.fai  1  fe )"  ' 

which  can  be  treated  similarly.  The  joint  cumulant  generating  function, 

°Q  n 

,  q2  )  =  iog  <p{qx ,  ?2  )  =  Z  Z  "  'fai  J  fe  )”  ' 

where  are  the  joint  cumulants,  can  be  used  for  faster  convergence  and  /r, can  be  expressed 
through  the  joint  moments  /u,  n_,  =  (rjl^r]n~1^ 


u=  Mu=\rl(rx)rl(r2))  =  cj-C 
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where  C 


n  -  7^2 


stands  for  a  two-point  correlation  function. 
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Figure  3.  Probability  density  distributions  of  the  sea  surface  cumulants  k2  ,  k,  /  cr3,...,  kX6  /  cr16. 
The  cumulants  of  a  Gaussian  variable  have  a  value  of  0,  indicated  by  the  dashed  vertical  line. 

RESULTS 

The  cumulant  expansion  leads  to  the  following  form  of  the  sea  surface  scattereng  cross-section. 


The  first  factor  exp(-2crg2)  corresponds  to  the  case  of  strictly  Gaussian  ocean  surface,  as  considered 

by  Ament  (1953)  and  Eckart  (1953b).  The  factors  that  follow,  represent  correction  to  the  scattering 
cross-section  associated  with  the  surface’s  deviation  from  Gaussianity.  The  signal  amplitude  (and 
intensity)  are  determined  by  the  factors  corresponding  to  even-order  cumulants.  The  factors 
corresponding  to  odd-order  cumulants  affect  only  the  phase  of  the  scattered  signal.  Nevertheless,  these 
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odd-order  factors  are  important  for  determining  the  nulls,  i.e.  the  locations  where  the  scattered  signal 
interferes  destructively  with  the  signal  from  the  source  and  thus  produces  “blind  spots”. 

For  the  case  of  incoherent  scattering,  the  cumulant  expansion  predicts  the  following  intensity 


Kc  =  exp[-  a-q1  (l  -  C)]x  exp 


•+ 

““(*■4,0  -  4/^34  +  3at2  2 ) 


x  exp 


q 

360 


(^6,0- 6^5,1  +15^4.2-10^3,3) 


Again,  the  first  factor,  exp[-crg2(l-C)J,  corresponds  to  the  case  of  strictly  Gaussian  surface,  and  the 
factors  that  follow  are  corrections  needed  to  account  for  the  deviations  form  Gaussianity. 

IMPACT/APPLICATIONS 


The  corrections  for  non-Gaussian  sea  surface  statistics  become  more  pronounced  with  the  increase  of 
the  parameter  [cr^q).  For  small  ( <j  q )  the  cumulant  analysis  predicts  changes  in  the  phase  of  the 

scattered  signal,  while  changes  in  the  intensity  remain  negligible.  Nevertheless,  it  is  reasonable  to 
assume  that  a  radar  signal  propagating  in  strong  refractive  duct  will  encounter  increased  (as  compared 
to  a  weaker  duct)  number  of  scatterings  from  the  surface  and  even  small  non-Gaussian  corrections 
could  accumulate  (over  multiple  acts  of  scattering)  to  a  substantial  change  in  the  signal’s  intensity. 

I  am  working  in  close  collaboration  with  Amalia  Barrios  of  SPA  WAR  to  explore  the  effect  of  non- 
Gaussian  sea  surface  statistics  on  the  predictions  from  the  Advanced  Propagation  Model  (APM).  Very 
preliminary  results  indicate  that  for  S-band  signals  (i.e.  longer  wavelengths  and  [crtlq)«  1)  the 

cumulant  corrections  are  insignificant,  for  intermediate  wavelengths  (X-band  signals)  the  cumulant 
corrections  visibly  improve  the  agreement  between  model  predictions  and  direct  observations  (Figure 
4).  For  shorter  wavelengths  (Ku-band  signals)  the  discrepancy  between  model  and  observations  is  also 
reduced,  but  with  some  increased  standard  deviation  of  the  difference  between  them. 
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Figure  4.  Comparison  of  direct  observations  from  RED  experiment  with  predictions  from  the 
Advanced  Propagation  Model  with  different  scattering  cross-sections  (Anient  (1953),  Miller,  Brown, 
Vegh  (1984),  and  cumulant  analysis)  for  X-band  signals.  Calculations  were  performed  and  the 
figure  was  provided  by  Amalia  Barrios  of  SPA  WAR. 


TRANSITIONS 

A  logical  transition  from  this  work  is  an  experimental  testing  of  the  available  scattering  models.  Data 
from  direct  measurements  suitable  for  such  testing  so  far  are  scarce.  Karasawa  and  Shiokawa  (1988) 
have  reported  measurements  intended  to  examine  the  predictions  of  Miller,  Brown, Vegh  (1984),  but 
do  not  reach  a  definitive  conclusion  (Figure  4  in  Karasawa  and  Shiokawa  (1988)). 
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